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We show how the Landau-Pomeranchuk-Migdal effect on photon production rates in a quark- 

Qh! gluon plasma can be derived via the usual Boltzmann equation. To do this, we first derive the 

electromagnetic polarization tensor using linear response theory, and then formulate the Boltzmann 

equation including the collisions mediated by soft gluon exchanges. We then identify the resulting 



X 

C^ ' expression for the production rate with that obtained by the field-theoretic formalism recently 

proposed by Arnold, Moore and Yaffe. To illustrate the LPM effect we solve the Boltzmann equation 
in the diffusion approximation. 
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The Landau-Pomeranchuk-Migdal (LPM) effect [1,2] plays an important role in the 
diagnostic tools of quark-gluon plasmas formed in ultrarelativistic nucleus-nucleus colli- 
sions, both in the energy loss [3], and the emission of photons and dileptons from the 
plasma, our focus here (for a recent review see [4,5]). The LPM effect takes into account 
multiple scatterings of the emitters and subsequent interference of the emitted radiation 
leading to a suppression of the bremsstrahlung rate from that obtained by Bethe and 
Heitler [6]. 

In recent discussions [7], the LPM effect has been calculated by explicitly summing 
the infinite series of Feynman diagrams that correspond to multiple scatterings. While 
providing a modern field-theoretical derivation of Migdal's results, such an approach is 
complicated by the need first to identify the relevant series of diagrams, and then to 
approximate these diagrams carefully in order to obtain useful expressions. As in the 
Landau theory of Fermi liquids [8] , multiple scattering processes involving a sequence of 
singular denominators are most effectively dealt with in the framework of the Boltzmann 
equation: not only does the Boltzmann equation capture the relevant diagrams, it also 
has the necessary kinematical approximations for small-angle scatterings built in via the 
gradient expansion of the collision term. However, neither the field theoretic treatment 
of [7], nor in fact Migdal's original derivation, make manifest the fact that the effect of 
multiple collisions is entirely captured in the relevant kinematical regime by the usual 
linearized Boltzmann equation. ^ Establishing this simple connection is the main purpose 



^ The first and third references in [7] do note an integral equation analogous in structure to a 
Boltzmann equation, but this analogy is not explored further. The Boltzmann equation derived 
here should not be confused with that discussed in [10]: in the latter paper, the LPM effect enters 
as a correction to the 1^2 and 2^1 collision terms. Migdal's original derivation describes the 
propagation of a charged particle interacting with fixed scattering centers at random locations. 
In a sense, the derivation presented in this paper extends Migdal's work [9] to the case where 



of the present Letter; we derive the LPM effect based on a (hnearized) Boltzmann equa- 
tion. While the rates obtained by solving the Boltzmann equation are not different from 
those obtained by the formalism of [7], the use of the Boltzmann equation significantly 
simplifies the derivation of the LPM effect and provides a new, more intuitive, perspective 
on the problem. Moreover, it is also the basis of new tools for addressing this issue in a 
non-equilibrated plasma. 

To leading order in the electromagnetic fine structure constant, a, the photon pro- 
duction rate is [11, 12]: 

Here u = |k|, and njjj,(c^,k) is the Fourier transform of the finite temperature current- 
current correlation function {K = (c<j,k)): 

n<,(c., k) = / d^X e^^-^ {U0)UX)) , (2) 

with j^{X) = etlj{X)j^ip{X) the electromagnetic current and X = (t, x) denotes the 
space-time coordinates. We use a metric with g'^^ = 1. To derive (1) we use the transver- 
sality of n^j^(-ft') : k'^Il1^^{K) = 0. A similar formula exists for lepton pairs for which 
K^ = o;^ — k^ > 0. For real photons, one can replace g^'^U^^ by g^^Hfj, where i,j = 1,2 
are the two directions transverse to the photon momentum; indeed the transversality of 
n^^ ensures that non-transverse polarizations do not contribute. The brackets in Eq. (2) 
denote a thermal average. The correlation function n'^j^(cj,k) is related to the retarded 
electromagnetic polarization tensor through (see, e.g., [13, 14]), 

n<,(c.,k) = --7I— imn-*(^,k) . (3) 



scatterings are due to two-body collisions. The present generalization focusses on soft photons, 
while Migdal's approach is valid for arbitrary photon energies. 
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Fig. 1. Left: quark-quark elastic scattering in the Born approximation via gluon exchange. The 

straight hnes denote quarks, and the curly lines gluons. Dynamical screening corrections, denoted 

by the thick dot, can be included via the hard thermal loop expansion. Right: the cut on the 

internal quark loop bubble that is used in connection with Eq. (17) below. Note that the scattering 

partner, as well the excitation running in the cut loop, can also be a gluon. 

For uj <^T the photon production rate per unit volume and frequency is thus 

^--|.i:-n.'fKk). (4) 

The main task in estimating the photon production rate is therefore to calculate nj|*(c(j, k). 
Since we can write the self-energy in terms of the response of the electromagnetic current 
to an external electromagnetic vector potential, 

the problem reduces to calculating (j*) in the presence of an external Aj. 

As we show below, the physics of the LPM effect is included in a calculation of the 
response of the current to an external field via the Boltzmann equation, even with the sim- 
ple collision term describing scatterings at the Born approximation level (see Fig. 1). The 
solution of the Boltzmann equation takes into account repeated scatterings - processes 
needed to include the physics of the LPM effect - as illustrated in Fig. 2. 

The leading order contributions to the photon rate, of order ag (= g'^/An, with g the 
strong coupling constant), correspond to real gluon-photon Compton scattering {qg — > 




Fig. 2. Resummation of ladder diagrams in the photon polarization tensor taken into account 

by the Boltzniann equation. The thick dot denotes dynamical screening corrections. 

g7 or qg -^ q'-f) and quark-antiquark annihilation {qq -^ g'j). These processes can be 

calculated by including hard thermal loop (HTL) [15, 16] corrections in the propagators 

[17, 18], and do not require further resummation. In particular they are not affected by 

the LPM effect. Since we focus here on the LPM effect we omit out these processes in the 

following discussion. 

The processes shown in Fig. 3, formally of next order in as, are collisions involving 
space-like gluons. The one-loop correction in the gluon propagator shown in these three- 
loop processes is only the first correction; the full correction should be carried out in 
practice by an HTL resummation. Naive power counting suggests that these two diagrams 
contribute in 0{a'^); however, "collinear enhancement" turns the contribution of these 
diagrams into a contribution of order ag- 
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Fig. 3. Higher order processes that are promoted to 0{as) by collinear singularities. 

We illustrate the origin of this enhancement by first studying the extent to which the 
quark of momentum P + K, between the quark-gluon vertex and the photon emission 
vertex in Fig. 3, is off-shell. On-shell, P^ = rn^, where rrioo ~ gT is the thermal mass of a 
quark of momentum ~ T. To estimate the virtuality of the intermediate quark, we work 
in the frame in which the photon four-momentum is K = [uj, 0, 0, kz = oo)', then 



{P + KY-ml, = 2P-K = 2u{^Jpl + ml-p,) (6) 

where m\ = p^ + m^. The right side of Eq. (6) becomes very small when m\ <^ p^, 
as occurs for small mass and emission of the photon in the forward direction, p_L ^ 
(collinearity). In this limit, the diagrams in Fig. 3 become singular. The quark thermal 
mass, which arises from HTL resummations on the quark lines (not explicitly shown in 
Fig. 3), prevents these diagrams from being truly singular, but the region of phase space 
where the quark and the photon are nearly collinear leads to a contribution ~ T^/m^ ~ 
l/cts. Combining this contribution with the explicit a^ from the vertices, we see that these 
diagrams become 0{as). 

K 
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Fig. 4. A virtual quark of momentum P + K emitting a real photon [K"^ = 0) and an on-shell 
quark of momentum P (P^ = m^ '' 



"oo/ 



In fact, a similar collinear enhancement affects an infinite set of processes. The 
collinear enhancement in the diagrams of Fig. 3, due to the small virtuality of the quark 
that emits the photon, can be rephrased physically in terms of a large photon formation 
time, or equivalently, the small energy denominators in the intermediate states. For the 
process in Fig. 4, the formation time is tp = ^/5E, with 



6E = uj + e^-e^+i,^— ■ — -, 7) 

2 Pz{Pz+^) 



where e^ = Jml^ + p^, and we assume m_L <^ p^- Typically, in a quark- gluon plasma, 
m^ ~ gT, while Pz ~ T. Thus for a photon of energy cj ~ T, we have 5E ~ g^T. But g'^T 
is in fact the characteristic scale of the rate of collisions with small (~ gT) momentum 
transfer Q. To see this result we write the scattering cross section as a = / dQ"^ {da / dQ"^) , 
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where typically da/dQ"^ ~ Q^/Q'^- Thus the collision rate, 7 = na, is ~ g'^T^ JdQ^/Q"^, 
where we use n ~ T^. As we verify later (see the discussion after Eq. (22)), the LPM effect 
regulates the integral in such a way that it becomes infrared finite but remains dominated 
by the contributions of soft momentum transfers Q ~ gT. It is thus of order l/{gT)'^, 
leading to the finite result 7 ~ g'^T^ /{gT^ ~ g^T. The formation time of a photon of 
energy a; ~ T is thus of the same order of magnitude as, or larger than the quark mean 
free path between two soft collisions, i.e., tp ~ I/7, the collision time. The formation time 
of soft photons uo ~ gT is even larger. 

Under such conditions, effects of multiple collisions on the production process cannot 
be ignored. Multiple scattering reduces the rate compared to that were all collisions treated 
as independent sources of photon production - the LPM effect. The multiple scattering 
diagrams that must be resummed in the polarization tensor are the ladders in Fig. 2. 
These processes, together with the self-energy corrections that need to be included on the 
quark lines, are the typical diagrams taken into account by the Boltzmann equation [13] 
(see also [23]; for a recent derivation in the context of QCD, see [24]). 

We turn then to the explicit formulation of the photon production rate using the 
linearized Boltzmann equation. The state of the system is described in terms of the dis- 
tribution functions of charged particles, which we denote by n/ for quarks and nj for 
antiquarks of flavor / (in order to simplify the discussion we ignore the gluons, on which 
the charged particle can scatter; including their contribution poses no conceptual prob- 
lem). Our task is then to determine the n's for a system initially in equilibrium perturbed 
by a weak electromagnetic potential, Aj\ then Uf = n^r + 6nf where n'r is the equilibrium 
distribution function for quarks of flavor /; similarly, Uf = fA + 6nf = n? + 6nf. We 
calculate the (5n's explicitly from the linearized Boltzmann equation. For soft photons, 
ko,k <tiT, this equation takes the form. 



{v ■ dx)Snf{p, X) + Cf^ ■ E— = C[5nf, 5%]; (8) 

here v^ = (1, v), with v = p/cp, so that v ■ dx = (9^ + v ■ V. In the force term, E = 
— VA° — dA/dt is the electric field, and ej is the charge of a quark of fiavor /. (For 
initially isotropic distributions, the magnetic field does not contribute to the force acting 
on the particles, to lowest order.) The 6nf are governed by a similar equation with ej 
replaced by — e/. The collision term, C, on the right side of Eq. (8) is linear in the 5n's. 
The electromagnetic current (j*) is given in terms of the (5n's by 



if) (X) = 2iV, Y: ^fj 7^^^ (^Mp, X) - %(p, X)) , (9) 



where Xc = 3 is the number of colors, and the factor 2 accounts for the two spin states of 
the quarks. 

The solution of the linearized kinetic equation is proportional to the forcing term ~ E 
on the left. It is convenient to write the (5n's in the form, 

5n;(p,X)^-e^W-(p,X)^ = ^W-(p,X)4(ep)(l-n;(ep)), (10) 

with the same equation for (5n/(p, X) with e/ — > — e/. The deviation W, the same for 
quarks and antiquarks, can be interpreted in terms of a distortion of the local momentum 
distribution caused by the shift de-p = — e/W^(p,X) of the single particle energies [19]: 
n/(p,X) = nj(ep) + 5n/(p, X) = nj(ep — e/W^). In terms of W^(p,X), the kinetic equation 
reads 

vdxW{p,X)-vE = C'[W], (11) 

with 

C = -ef{dnydep)C'. (12) 



Fourier transforming with respect to the spatial coordinates we rewrite the kinetic 
equation as 

ivKW{p,K)+v-E{K) = -C'[W], (13) 

where K is the four-momentum of the produced photon, and we use the same symbol for 
a function and its Fourier transform, e.g., W{p,X) and W{p,K). Using Eq. (10) in (9) 
we find, after Fourier transforming, 

_2 f d^P i,,,, j.Jn} 



where e^ = AN^ J2f ^}- The extra factor 2 in Eq. (14), as compared to (9), accounts for the 
equal contributions of quarks and antiquarks. We calculate the polarization tensor from 
Eqs. (14) and (5) below. 

The linearized collision term, with Eq. (10), reads 



|2 



C = -%Y.I (27r)V^)(P + Pi -P,- P3) I^P-P^-P3.P^i 

xn;(p)(l - n;(p3))n°,(pi)(l - n],{p^)) [W{p, K) - W^(p3, K)] . (15) 

Here all quarks are on their mass shells, and /p. = /d'^pj/(27r)^. The matrix element 
squared, |A^p,pi^p3,p2 1 , is that for one-gluon exchange, as depicted in Fig. 1; it is averaged 
over the spin and color states of the incoming particle (of momentum P) and summed 
over the spin and color states of the other particles. To obtain (15) we also use the fact 
that the terms involving 5nfi{pi^K) and 5nfi{p2-,K) cancel when summed over quarks 
and antiquarks (e.g., 5nfi{Tpi,K) + 5nfi(j>i,K) = 0). 

At this point, it is convenient to use the four- momentum transfer Q = P2 — Pi as an 
integration variable. One can then perform the integrations over p2 and ps, and obtain 
the scaled colhsion term (12) as 
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C' = - E / (^27r5(go - V . q) [W^(p, K) - W{v - q, K)] 
X / 2^5{q, - vi ■ q)4(pi)(l - 4(Pi)) '^'''^^'"'''^^''' ' (16) 



where v = p/cp and Vi = pi/ep^, and we have used the fact that the momentum transfer 
is small (g <^ pi) in order to simplify the delta functions and the statistical factors, e.g., 
writing ep_|_q — ep ^ v ■ q, and n'r(p — q) ~ n^A'p). After these simplifications, the integral 
over pi can be done and the contribution of the matrix element expressed in terms of the 
spectral density p^^, {Q) of the gluon propagator in the HTL approximation (see Fig. 1, 
and Ref. [20] for details). One then arrives at the scaled collision term 

C'= -g'Cj J ^,S{qo - q ■ v)^^v'^v^p;:\Q)[W{p, K) - W{p - q, K)], (17) 

where C/ = (A^^ _ i)/2N^. 

In order to proceed further, we examine the specific angular dependence of the fluctu- 
ation W^(p, K) involved in the emission of soft real photons (typically with u ~ gT ^ T). 
We first note that by symmetry the solution of Eq. (13) with (17) must be of the form 

l^(p,ir)=vE/(p-k,p), (18) 

where E is the unit vector along the direction of the electric field. This angular structure 
is illustrated by the collisionless Boltzmann equation, (13) with C = 0, which has the 
solution 

mw X v-E(K) , , 

IV ■ K 

For real photons, f-A' = c<j(l — v-k), with k = 'k/uo. In the case of massless particles, 
V ■ K vanishes when p is parallel to k, leading to a diverging iy'^°)(p, K) - the "collinear 
enhancement" discussed earlier. Indeed, for soft photons the drift term is simply the 
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energy difference 5E defined in Eq. (7), i.e., v ■ K k, oom\/2p'^ (with m?^ = p5_ + ^ooi 
p_L _L k, and p ■ k ^ p_L, moo)- Thus, noting that E(K) = iujA{K), we have: 

W.(o)(p K) = A^^(^^o^<P ^cos^ 

^^' ' 2(l-cos^)+m^/p2 - ^2 + ^2^/p2' y") 

for small ^, where 6 is the angle between p and k, and the angle between E and p. In 
the case of massless particles, Eq. (20) exhibits the small angle divergence W^'^^ ~ 1/^ 
mentioned above. For massive particles, W^'^^ vanishes at ^ = 0, but remains peaked at 
small 6 ^ 00 = rrioo/p ~ g- 

Collisions, dominated by small angle scattering, maintain the peaking of the LPM 
fluctuations Vr(p, K) at small forward angles. Accordingly the solutions Vr(p, K) of the 
linearized Boltzmann equation are of the form W{p,K) = v ■ Fih{6,p), with h strongly 
peaked at small 6. This structure simplifies the calculation of the collision term, as we now 
show, and confirms the kinematical approximations that we made in deriving Eq. (17). 
Since in a collision |p — q| differs from |p| by subleading terms, the magnitude of p re- 
mains basically constant during collisions, with the direction of the dominant p remaining 
approximately aligned with the momentum of the photon. Because q is primarily trans- 
verse to p, we can neglect the dependence of W on q^. Thus q ■ p ~ Q'z, which allows us 
to integrate over go and qz in Eq. (16). Using the sum rule in [22] we find, 

f dqodqz v'^v" htl 1 1 .^, ^ 

J 2tt go q± q± + "^^ 

from which we obtain, 



C = -g'CjT I 7T^^7^^^3:7t[W^(p) - W(p - q^)] . (22) 



(fqj_ m^ 

{27ryql{ql + ml] 

We do not explicitly indicate the dependence of W on K. Since W^(p) — W^(p — q) vanishes 
smoothly as q ^ 0, the integral in Eq. (22) is infrared convergent. We expect W^(p — q) 
to decrease rapidly when \p± — g±| ^ moo, as in the coUisionless case, Eq. (20). Thus the 
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integrand in Eq. (22) is dominated by momenta of order rriD ~ ttIqo ~ qT ^ and expect the 
integral to remain of order g^T. ^ 

We now show that the present Bohzmann equation, with the colhsion term (22), leads 
to the formulation of the LPM effect of Arnold et al. [7]. We first observe that the kinetic 
equation that appears in [7] is for a vertex function rather than for a particle distribution. 
The quantity f * of [7] is related to the present W by 

Jam = ~%^ ^^"-^ ' ^ ^ 

where * denotes the complex conjugate. In terms of /, 

.mn;i.W)^|lf*,MRe/|EL.|.;. p4) 

which coincides with Eq. (2.1) of the second of Refs. [7]. Furthermore, the function f* 
obeys the kinetic equation we obtain by taking the functional derivative of the linearized 
Boltzmann equation with respect to Aj. This equation is identical to Eq. (2.2) of the 
second of Refs. [7] , after identification of the energy difference 5E with the drift term in 



^ Although the electrical conductivity can be obtained as the limit of the polarization tensor as 
w, k ^ 0, we cannot directly use Eq. (22) to derive this limit. The peaking at small angles in the 
LPM effect is in contrast to that in transport calculations [21], such as the electrical conductivity 
[25, 26]. There, a similar cancellation of the small q^ contributions in the term ll^(p) — W{p — 
q) makes the collision integral converge in the infrared, but the angular dependence of the 
fluctuation VF(p) induced by a uniform electric field does not constrain the momentum transfer 
to be soft, and we may not neglect the qz dependence of W as in deriving Eq. (22). The collision 
term involved in the conductivity is proportional to g'^T ln{T / gT) , where the upper cutoff in the 
logarithm comes from the limit of validity of the soft momentum approximation (q ^T), while 
the lower one originates from screening [25]. We defer discussion of the transport solutions of 
the Boltzmann equation to a future publication. 
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our Boltzmann equation (see the discussion below Eq. (19)). 

As an illustration of how the LPM effect emerges from the present Boltzmann equa- 
tion, we write, following Migdal [2], a diffusion approximation for the collision term (22). 
This approximation is not quantitatively useful for describing quark-gluon plasmas, since, 
as we shall see shortly, it requires rrioo ^ itlo, a condition which is not realized in quark- 
gluon plasmas where rather rrioo ~ ""^d (^^L — fi'^^^C*//4, and mj^ = {2Nc + Nf)g'^T'^/6). 
However, it yields simple analytical expressions which allow us to illustrate certain of the 
physical points made in the preceding discussion. The diffusion approximation is derived 
by expanding W{p — q^±) to second order in q^: 

where the diffusion constant is D = {g'^CfTm'jj/SiTp'^) \n{qmax/'m'o)y &^d Qmax ~ ""^oo- That 
rriao is the appropriate upper cutoff can be understood from the following argument. Since 
H^(p) decreases rapidly when p± ^ rrioo (cf. Eq. (20)), we conclude that the integrand 
in C in Eq. (22) is approximately constant when qj_ ^ rriD and it behaves as l/q'j_ for 
fn^ <^ q_L <^ rrioQ, and as l/g| for qj_ ^ moo. Thus qj_ ~ moo is the appropriate upper 
bound for the integration over q±, and the expansion of the collision term involved in 
the diffusion approximation may be viewed as an expansion in m^/moo- Note that the 
assumption m^ <^ moo justifies the leading log approximation used in estimating D. 

The Boltzmann equation in the diffusion approximation then reads 



y (vx' + ml/p')Wip) + DVl^Wip) = -luw ■ A. (26) 

Again, by symmetry the solution is of the form ^^(px) = v_|_ ■ Ay9(s)/s, where s = 
v^/2 ~ (sin^6')/2. The equation for ip{s) is 
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d 

iuj{s + SQ)ip{s) + 2sD -—ip{s) = —ius, (27) 

with So = '"^L/2p^- In terms of (/?, the current (14) is 



_2 / <i^P dn'^ / \ At ^^ /"°^„2j„'^^° 



where we use k ■ j = k ■ A = 0, and the fact that ip{s) is a rapidly decreasing function of 
s, in order to extend the range of the s-integration to +oo. From Eqs. (28) and (5) (with 
gij = —Sij) we derive 

e^ f°° driQ 



27r^ Jo dp 

where ^{uj;p) = J^ dslunp{s); the dependence oi ip on p (not indicated exphcitly) comes 
from the diffusion constant D ~ 1/p^ and Sq = m^/2p^ in Eq. (27). 

Solving Eq. (27) by iteration, we derive the solution as an expansion in powers of D, 
i.e., in the number of collisions. In zeroth order (no collisions), v?*-"-' = —s/{s + Sq), which 
is real and does not contribute to $: there is no radiation in the absence of collisions. 
Substituting the lowest order result into the collision term, we find the single collision 
contribution, ip'^'^^ = AiDssQ/(uj{s + sq)^)) which is imaginary and yields $*-^-' = 2D/3ujSo, 
independent of p. Note the role of the fermion mass, entering through the factor sq; as 
^oo -^ 0, $*-^-' diverges - the coUinear divergence discussed earlier. Using this expression 
for $*^^) in Eq. (29), we recover the low frequency Bethe-Heitler rate from Eq. (4): 



'''?•' - ^^>,. ^ C./A C. - ^^'4 in ^. (30) 



d^x 127r2 ^" uj ' " 727r3 m 



oo 



It is easily verified that this expression agrees with that obtained from the formulae in 
Sec. 4.2 of Ref. [22] (in the appropriate limit, vrtao ^ ^^d\ an a postenon justification for 



our use of qmax = i^oo in Eq. (25). The emission rate is of order g^ and falls as \/m 



2 

oo' 
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Proceeding further, we find that $'^^) = and $*^^^ ~ D^/(uj^sl). The iterative solution 
breaks down when $*^^^ ~ $*^^\ which occurs when to ^ D/sq. In the small u regime, we 
must solve Eq. (27) more accurately. The exact solution is in fact known; since Eq. (27) is 
identical to Eq. (44) of Ref. [2], we exploit the analysis there write the exact solution in the 



form $(c<j) = {2D/3soUj)(f){T), where r = {sq/A)Juj/ D. The function 0(r) satisfies 0(r -^ 
0) ^ 6r and 0(r ^ cx)) = 1. Thus, when u > D/sl = D2p^/ml^, $(cj) ^ 2D/?,Squj, and 



one recovers the Bethe-Heitler limit. On the other hand, as c<j ^ 0, $(c<j) ~ J D/uj, and 
we see that the rate is suppressed by a factor ~ ^/cJ, the LPM effect. Then 

dNi^^'' [ZidN^^ du 

a^x y ujc d^x v^jcjc 

where 

162(ln2)2m^ \m^ J ^ ^ 

Here we obtain the particular form of the spectrum ~ dioj^/Zj from the diffusion approx- 
imation. However, the same form of the spectrum emerges, as we find, from a numerical 
solution of the Boltzmann equation in the regime, rrioo ~ ^^d, where the diffusion ap- 
proximation is no longer valid. The diffusion approximation primarily affects the overall 
normalization, and has little effect on the shape of the spectrum for small photon energies. 

In this Letter we have considered only the case of soft photons. More generally - and 
in particular for hard photons with cu ~ T - neither the energy difference 6E^ nor the 
coupling between the quarks and the applied electromagnetic field, can be approximated 
in a gradient expansion. However the kinematical conditions that allow one to obtain 
the linearized collision term from the quantum field equations still hold. The resulting 
Boltzmann equation takes a similar form, but with more accurate drift and Vlasov terms 
(see, e.g., [13]). 
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